Antiproton-Proton Channels in J/ip Decays 
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The recent measurements by the BES Collaboration of J/ip decays into ypp indicate a strong 
enhancement at pp threshold not observed in the decays into -K°pp. Is this enhancement due to a 
pp quasi-bound state or a baryonium? A natural explanation follows from a traditional model of 
pp interactions based on G-parity transformation. The observed pp structure is due to a strong 
attraction in the 1 So state, and possibly to a near-threshold quasi-bound state in the 11 So wave. 
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A point of interest in the antiproton interactions is the 
question of existence or non-existence of exotics in the 
nucleon-antinucleon (NN) systems: quasi-bound, vir- 
tual, resonant, multiquark or baryonium states 0. Such 
states, if located close to the threshold, may be indi- 
cated by large scattering lengths for a given spin and 
isospin state. For this purpose the scattering experiments 
are apparently the easiest to perform with a good pre- 
cision. However, a clear separation of quantum states 
is not easy. Complementary measurements of the X ray 
transitions in antiprotonic hydrogen are useful to select 
some partial waves. These are particularly valuable when 
the fine structure of levels is resolved. Such a resolution 
has been achieved for the IS states Q and partly for the 
2P states Another method to reach selected states 
are formation experiments. In this way a resonant-like 
behavior was recently observed by BES Collaboration in 
the radiative decay 



J/ip — > jpp 



(1) 



close to the pp threshold Q . On the other hand, a clear 
threshold suppression is seen in the pionic decay channel 



J/i/> 



n°pp. 



(2) 



To understand better the nature of the enhancement, 
one should look into the pp sub-threshold-energy region. 
This may be achieved, indirectly, in the pd low-energy 
scattering or pd atoms. Such atomic experiments were 
performed, but the fine structure resolution has not been 
achieved so far p|. 

The purpose of this letter is to discuss the physics of 
slow pp pairs produced in the J/ip decays. The J PC con- 
servation reduces the number of pp final states to several 
partial waves. These, denoted by 21+1 2S+1 Lj, differ by 
their isospin I, spin S, angular momenta L and total 
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spin J. Close to the pp threshold, quite different be- 
havior of scattering amplitudes is expected in different 
states. In the IS state of antiprotonic hydrogen it is 
the = ( n 5 + 31 S )/2 and 3 Si = ( 13 5i + 33 Si)/2 
waves which are studied |2(. While atomic experiments 
determine the scattering lengths, the BES experiment al- 
lows to extend this knowledge into a broad energy region 
above the threshold. As will be shown, the radiative J/ip 
decay involves also the 11 So + 31 So combination. The 
understanding of this and other involved states should 
be based on the experience gained in studies of elastic 
and inelastic NN scattering. We use the Paris potential 
model [E B f° r this purpose. 

To o ur p resent knowledge, none of the available related 
works [lfl [H[ll EUll on the J/ip decays (JTJ) and © 
have given a comprehensive explanation of the BES ex- 
perimental spectra. Only two of these papers [ToL IT^ 
compare their results to the data. The Jiilich NN model 
is used in 0] to show that, within the Watson-Migdal 
approach, the isospin 1 S-wave can reproduce the low 
energy part of the pp spectrum in the radiative decay. 
The same spectrum is fitted with a constant scattering 
length in Ref. The length obtained in this way is 

larger than the lengths calculated in potential models. 
In Ref. 1 more realistic but spin averaged constant 
lengths are shown to generate some low-energy enhance- 
ment in reaction ffl. A X-matrix, calculated with the 
one-pion exchange in the Born approximation, is consid- 
ered in Ref. . An enhancement is seen, but this model 
is too simple to describe the NN interactions. The for- 
mation mechanisms in the radiative decays are discussed 
qualitatively in Ref. fl3j | , where the quantum numbers of 
final states are listed with the recommendation to look 
into decay modes of the pp systems. 

In the present work the following results are obtained. 
The set of allowed final pp states is limited to three par- 
tial waves in the photon channel and to two waves in the 
pion channel. Among the three possible pp states in the 
ppy channel one is dominated, at very low energies, by 
the well known pp( 13 Pq) resonance, formed as a result of 
attractive one-pion exchange forces. However, this state 
as well as another allowed 3 Pi state cannot explain the 
experimental spectrum. The final ppj state is dominated 
by the pp^So) wave. A strong attraction arises in this 
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TABLE I: States of low energy pp pairs allowed in the J/ip — > 

7PP and J/t/j — > 7r°pp decays. The first column gives the 

decay modes to the specified internal state of the pp pair. 

Well established, two particle analogs are indicated in the 

second column jl5|. The third column gives J PC for the light 

spectator particles, photons or pions. The fourth column gives 

J PC for the internal pp system; the last column gives the 

relative angular momentum of the light particle vs. the pair. 
JPC = 1 ( — ) for jy^_ 



Decay mode 


Analog 


J PC [~f or 7T ] 


J pc [pp] 


Relative 1 


JPPCSo) 


777(1444) 


I— 


0-+ 


1 


-YPP( 3 Po) 


7/0 (1710) 


I— 


0++ 





~1PP{ 3 Pi) 


7/1 (1285) 


I— 


1++ 





•K°pp{ 31 Pi) 




0-+ 


1+- 





n°pp( 33 Si) 




7T p 


o-+ 


1~~ 


1 



wave as a result of coherent one-pion and two-pion ex- 
change forces. It produces broad, deeply bound states, 
difficult to detect. However, the recent version of the 
model H, adapted to hydrogen atom data, generates a 
near-threshold state in the related pp^So) wave. This 
state is about 50 MeV wide and bound by 5 MeV. 

In the ppir decay channel, isospin being conserved, 
two pp waves, 33 Si and 31 Pi, are allowed. These indi- 
cate distinctly different threshold behavior. The S wave 
is ruled out by the experiment and the pp( 31 Pi) leads to 
a natural explanation of the BES spectrum. These find- 
ings can be unified in a qualitative model for both decay 
modes. 

The allowed final states. The J PC conservation limits 
the number of slow pp final states. The latter are under- 
stood as pp pairs with small M p p — 2m p where M p p is 
the pair invariant mass. The allowed states are listed in 
Tableland a few possibilities exist for each channel. 

The BES experiment provides an angular distribution 
for the photons. The specific situation of e~e + collision 
is that the projectiles are polarized perpendicular to the 
beam direction and the J/ip spin direction follows that 
of the beam. The angular distribution can be measured 
within a limited range. It offers a hint which in princi- 
ple permits to further reduce the number of allowed final 
states. For completeness, the angular distributions of the 
photon in the states of interest is calculated below. The 
simplest Lorentz invariant couplings which are also gauge 
invariant are given in TablelTTl The pp pair is described as 
a single scalar ( 3 Po), pseudoscalar ( 1 S'o) or pseudovector 
( 3 Pi) particle. This picture is expected to work for the 
pairs of small center of mass (CM) system energies. Next, 
the appropriate couplings are reduced to the J/ip CM 
system and the corresponding angular distributions are 
calculated. The couplings given in Table ln| follow from 
Hamiltonians that couple the electromagnetic field f ap to 
the J/ip field F ap = d a V p — d p V^ and to the correspond- 
ing scalar ip s , pseudoscalar Lp ps and pseudovector ip^ 3 
fields describing the low energy pp pairs. These gauge in- 
variant Hamiltonians are H s = (en s /2Mj) J f ap F ap (p s , 



HP S = (eK ps /4Mj) / fapF^PipP 3 : 16j and H pv = 
e J frrpVa ' f 1 ^ £ apa>3 , where K s ' ps are the magnetic mo- 
ments for the transitions. The vector nature of the fields, 
expressed by V 3 ^ 11 = 0, leads to the results given in 
Table [HI The first two angular distributions are well 
known |l7| . 8 denoting the angle between the 7 emis- 
sion and the beam direction, the angular distributions 
(cos 2 8 + 1) and sin 2 8 were tested against the data in 
Ref. These have indicated a preference for the first 
choice, i.e. radiative transitions to 3 Po or 1 So states. 
However, as can be seen in Table UTI a transition to 3 Pi 
state is not excluded. 

Final state interactions. Any multichannel system 
can be conveniently parameterized by a i^T-matrix which 
guarantees unitarity of the description. The transition 
amplitude from a channel i to a two-body channel / may 
be presented in the form 



T if = 



A 



if 



1 + iqAf f 



(3) 



where Aif is a transition length, Aff is the scattering 
length in the channel /, and q is the momentum in this 
channel (see e.g. ^U). Both lengths can be expressed 
in terms of energy dependent K matrix elements. The 
same formalism describes the scattering amplitude in the 
channel / as 



T Sf - 



A 



If 



1 + iqAff 



(4) 



In the process of interest the formation amplitude Ai / is 
unknown, but Aff is calculable in NN interaction mod- 
els constrained by other experiments. For slow pp pairs 
the final state interactions in the ir°pp and jpp systems 
are dominated by interactions in the pp sub-system. A 
formal manipulation of Eqs. © and (0} yields 



™ At 

A ff 



A 



A 



If 



~± h (5) 



which defines a quantity Cif = Aifq L j Af f . For 5-waves, 
the standard final state dominance assumption (Watson- 
Migdal) is equivalent to a weak energy dependence in 
Cif. This is usually true in a small energy range where 
the denominator in Eq. provides all the energy de- 
pendence. In the pp states such an approximation is 
correct for q up to about 0.5 fm -1 . It fails at higher 
momenta since Aff is energy dependent. On the other 
hand Aif stems from a short range cc annihilation pro- 
cess. The annihilation range is of the order of l/m c pl| 
and only a weak energy dependence is expected in Aif. 
We assume Aif = Aif(Q)/ [l + (r q) 2 ] with a range pa- 
rameter r Q well below 1 fm. For P-wave final states, 
the low-energy behavior gives Aif « Ahq, Aff fa A^q 2 

where the Ah are parameters and the A 1 ^ the scatter- 
ing volumes. The latter are energy dependent as a result 
of medium ranged ir exchange forces. This dependence 
is particularly strong in those waves which involve res- 
onances. The Watson approximation is not appropriate 
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TABLE II: Radiative couplings of low energy pp pairs to 7 J/ip and angular distributions of the final photon with respect 
to the J/ip spin direction, J denotes the J/ip and h the pp pair. The second line specifies the transformation property of pp 
pair. The third line gives the invariant couplings and the fourth line presents these in the J/ip center of mass. k^' J ' h are the 
four momenta and e^ J ' h the polarization four-vectors of the corresponding fields. The last line gives the angular distributions 
for the three different cases. The photon energy is denoted by u>, that of pp pair by Eh and the J/ip mass by Mj. For 
E h = 2m p : u/E h = 0.46. 



decay mode 

Hpp) 

coupling 
coupling in CM 
angular profile 



7 PPCSo) 
pseudoscalar 
e^"klelkU J p /Mj 
(e J Ae 7 )-k 7 

1 + cos 2 e 



7 PPi 3 Po) 

scalar 
- k J p ei)FJ p /Mj 



we. e ' 

1 1 „„„2, 



7Pp( 3 Pi) 
pseudovector 

afiapui 7 J h 

(e J A e<) ■ [k 7 (k h ■ e h )/E h - e h to] 
2sin 2 6» + [(1 +Lj/E h )/{1 -u/E h )](l + cos 2 
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FIG. 1: 7r°pp( 31 Pi) decay channel. Experimental data have FIG. 2: ir°pp( 33 Si) decays. Data as in Fig. 1. (a) Final state 



2 of Eq. @. Constant AL(0) and formation 
0.55 fm are chosen to obtain a good 



been extracted from Fig. 2(a) of Ref. (a) Final state factor 
q | Tff/q | 2 (Watson approximation). Constant Cj/ of Eq. JjJ 
is chosen to fit the low-energy part of the data. Four versions 
of the Paris potential model jy, 0, 13 are used. This ap- 
proximation fails for M P p — 2m p > 40 MeV (g > 1 fm -1 ). (b) 
The rate q \ T if 
range parameter r, 
fit to the data. All four potentials give equivalent fits, even- 
though a 118 MeV wide state bound by 15 MeV is generated 
in version [f| in the 31 Pi wave. 



there, and Eq. J2J must be used. The transition length 

is parameterized as Aif = A}f(0)q/ [l + (r q) 2 ] 2 . 

The advantage of K matrix formalism is clear in the 
analyzes of low-energy final state interactions since it iso- 
lates the kinematic singularity into a definite form given 



factor q \ T ff | 2 of Eq. ©. (b) Rate q \ T if | 2 of Eq. ©. This 
wave is not consistent with the BES data whatever the choice 
of pp potential version, of Cif in (a) or Aif(0) and r in (b). 



by Eq. ©. The two functions Aff and Aif depend only 
on q 2 . Hence, close to the threshold, a constant scat- 
tering length approximation in Eq. JSJ) may well indi- 
cate some sub-threshold phenomena. This approxima- 
tion has been used in Refs. [Til IT^ . In the pp system 
the energy dependence in Aff is strong as pointed out in 
Ref. ^3 011 the basis of an one-boson exchange version 
of Bonn potential. A similar behavior is seen with the 
Paris model although these two potentials differ strongly 
in the two-pion sector. As shown below, Eq. JSJ with a 
constant C,/ and realistic Aff{q 2 ) describes a too narrow 
energy range. The selection of the best pp partial wave 
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FIG. 3: Jpp( S Pi) decays. Experimental data have been 
extracted from Fig. 3(a) of Ref. 4]. (a) Final state factor 
q | Tjf/q | 2 of Eq. JSJ. The enhancement in Paris 05 so- 
lution i s related to a 18 MeV wide state bound by 5 MeV 
in the 33 Pi wave, (b) Rate q \ T if | 2 of Eq. ©. This wave 
cannot reproduce the BES data whatever the choice of pp 
potential version, of df in (a) or A}f(0) and r Q in (b). 
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FIG. 4: 7Pp( 3 Po) decays. Data as in Fig. 3. (a) Final 
state factor q \ Tff/q | 2 of Eq. ©. The low energy part is 
dominated by the resonance in the 13 Po wave at 1876 MeV, of 
10 MeV width, present in all models. However, for q > 1 fm - 
this approximation fails to fit the data, (b) Rate q \ Tif | 2 of 
Eq. with r — 0.55 fm. This rate can describe only the 
q > 1 fm" 1 part of the spectrum. This wave is not consistent 
with the BES data. 



requires Eq. J2J|. This equation offers also an explicit and 
unique dependence on the the on-shell A tf (or Tf f since 
1/(1 + iqAff) = 1 - iqTff). An off-shell A ff may be 
involved in Ai f if one attempts to construct a model for 
the pp formation. 

Results. There exists substantial phenomenological 
control over Aff. Here these scattering lengths are 
calculated in terms of the Paris NN potential model 
and the same procedure is applied to both decay modes 
J — > Tr°pp and J/ip — > 7PP- Figures 1 to 5 present 
the results obtained with Eq. JSJ and Eq. (0 for the five 
states of interest calculated for the four versions of Paris 
model [a, 0, la, 0| which evolved over the last 20 years. 
This evolution followed the increasing data basis. The 
last version 6] is based on 3934 data which includes the 
recent antiproton-hydrogen widths and shifts 0, and 
the total np cross sections of Ref. [20J. The Coulomb 
interactions yield enhancements of the S waves at very 
low energies due to Gamov factors. These affect the final 
state interaction for q < 0.15 fm" 1 and produce spikes. 
Since the amplitudes are weighted by the phase space fac- 
tor q, these become unessential. The q factor represents 
a residual piece of the full three-body phase space ■ 

A qualitative model. As exemplified by the final state 



calculations, the BES findings are most consistent with 
a ppi^Pi) wave in the ir°pp channel and a pp( 1 5o) wave 
in the jpp channel. Therefore the experiment leads us to 
a simple picture for the slow pp formation. 

The initial heavy cc quarks in the J/ip state of 
J PC = 1 annihilate into a NN pair. As argued in 
Refs. [HEI that process is mediated by three gluon ex- 
change. Due to isospin conservation, the baryon pair is 
formed in an / = state of nfi+pp as indicated by exper- 
iment [2^ and calculations of Ref. • The pair inherits 
the J /if) quantum numbers J PC = 1 and forms a 3 Si 
state. Next, the emission of a pion or a photon takes 
place. The ttq emission proceeds via the standard irNN 
coupling (f^NN l^fn-ix) q ■ cr. It requires one nucleon to 
flip spin and change angular momentum, which leads to 
the final pp{ l P\) state. The photon may be produced as 
a magnetic one or as an electric one. The relevant for- 
mation amplitudes are given by the transition operator 
(e/2m p ) [2 e 1 ■ q + i er • (k 7 x e 7 )]. In the final states 
q is small. In the intermediate states it is not necessar- 
ily so, but any formation mechanism would favor small 
momenta. Since | k 7 | is large we conclude that it is the 
magnetic transition which is more likely to occur. It fa- 
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FIG. 5: jpp{ So) decays. Data as in Fig. 3. (a) Fi- 
nal state factor q | T// ] 2 of Eq. 0. At higher (q > 
2 fm - ) momenta this approximation begins to fail, (b) Rate 
q | T if I 2 of Eq. © with r a = 0.55 fm. The latest Paris 
model Q offers the best fit to the data with an 11 So wave 
involving a quasi-bound state located very close to threshold, 
of 53 MeV width and 5 MeV binding. 



vors formation of the final jpp^So) state which arises in a 
most natural way. In the initial pp( 3 Si) wave, the proton 
and antiproton magnetic moments are opposite and the 
transition to pp( Sq) involves spin and magnetic moment 
flips. Large moments create large radiative amplitudes. 
The emission model indicated above, yields comparable 
branching ratios of the 7 and 7r° channels, as found in 
the experiment. This ratio follows roughly the ratio of 



the coupling constants f^M N / (4 e 2 ) w 2.8 while the ex- 
perimental ratio is « 3 15]. 

The final pp state involves the isospin 1 plus isospin 
combination. The pair may be also formed in the nn state 
and undergo a transition to pp in the final state. That 
process is expected to be suppressed, since that transition 
implicates the Tff(I = 1) — Tff(I = 0) amplitude which 
is about an order of magnitude smaller than the elastic 



T ff (I = 1) + T f f(I = 0) one. The simple model of 
final photon radiation discussed above would reduce the 
neutron channel even further, due to different charges 
and magnetic moments. 

Conclusions. We have shown that the new results of 
the BES Collaboration find a natural explanation in a 
fairly traditional model of pp interactions based on G- 
parity transformation, dispersion theoretical treatment 
of two pion exchange and semi-phenomenological absorp- 
tive and short range potentials. This model predicts 
quasi-bound states close to the threshold, in particular 
in the pp( 33 Pi) and pp( 11 S'o) waves and a resonance in 
the pp( 13 P ) wave. The first two indicate a strong de- 
pendence on the model parameters and, so far, are not 
confirmed in other experiments. The third one, the res- 
onant state, is well established. 

It is the So state which reproduces the jpp spectrum 
found by the BES collaboration. This wave is domi- 
nated by a strong attraction due to the pion exchange 
forces. This attraction generates broad, deeply bound 
states. The recent atomic and scattering data indicate 
that such a state in the 11 So wave is located close to the 
threshold. The BES data offer some support for the ex- 
istence of such a state. The actual energy level and its 
width are affected by interactions at distances less then 
1 fm. These are not not fully understood and only partly 
controlled through phenomenology. 

In order to better see the nature of the 11 S'o state, one 
should look directly under the pp threshold. This could 
be done with measurements of the invariant mass of few 
meson systems coupled to pp just below the threshold. 
The selectivity in partial waves is necessary, and a con- 
venient way to reach that is the J/ip — > 7 mesons decay. 
Another, indirect method is to achieve a fine resolution 
of energy levels in antiprotonic atoms. Some anomalies 
were found in atoms with nuclei characterized by weakly 
bound valence protons j2^| • These anomalies may reflect 
a resonant behavior of the pp scattering amplitudes in 
the region of pp quasi-bound states. More systematic 
measurements are necessary to pinpoint the pp wave re- 
sponsible for these effects. 
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